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Heavy flavour production and spectroscopy are key components of the HERA physics programme. I will summarise 
a selection of the recent results obtained by the HI and ZEUS collaborations. The production of excited charm 
mesons and J/tp will be discussed as well as measurements of b quark cross-sections in photoproduction. The status 
of searches for exotic bound states and the D*p resonance will be updated. 



1. INTRODUCTION 

Heavy flavour production in e^p collisions at HERA provides a good testing ground of perturbative Quantum 
Chromodynamics (pQCD) as the high quark mass provides a hard scale. Furthermore, other hard scales such as 
Q 2 , the virtuality of the exchanged boson, or p t , the transverse momentum of the heavy quark, allow resummation 
techniques to be tested. Measurements of production rate and kinematic properties of heavy quark bound states 
such as charmonium also give direct access to the non-perturbative part of the production process. New theoretical 
models can also be tested by searches for exotic bound states like the D*p resonance. 

Different kinematic variables are used to describe the ep interaction at HERA: Q 2 , the Bjorken scaling variable, x, 
and the inelasticity, y. Until 1997 HERA ran at a centre-of-mass energy of y/s = 300 GeV. This energy was increased 
to y/s = 320 GeV for data taken from 1998 onwards. Due to improvements of the detector and accelerator during a 
break in the data taking, the available dataset is split into two periods. In the HERA I period from 1996 until 2000 
about 130pb _1 and between 2003 and 2007 of HERA II about 400 pb -1 per experiment were collected. 

The kinematic range of the analysed data can be separated in the following two regimes: photoproduction (7p), 
where the exchanged photon in the process is almost real, and deep inelastic scattering (DIS), where the exchanged 
photon is virtual. Experimentally, jp is defined by the scattered electron not being in the acceptance region of the 
detector, corresponding to a cut Q 2 < 1 GeV 2 . 



2. THEORY 



There are different approaches for the calculation of heavy flavour production in next-to-leading order perturbative 
QCD. The massive approach assumes no initial charm or beauty in the proton (or photon). Heavy flavours are only 
generated dynamically from the gluon distribution. This approach is particularly valid if the heavy quark mass is of 
the order of other hard scales like Q 2 or the transverse momentum p t of the heavy quarks. In the massless approach, 
the leading-logarithmic and next-to-leading logarithmic terms for example in a s logQ 2 /m 2 b are resummed. This 
approach assumes the heavy quark to be massless and is therefore only valid if other quantities like Q 2 or p t provide 
the dominant scale. Models using the kx factorisation approach [l| are based on non-collinear parton dynamics based 
on the CCFM Q evolution equations. In heavy quark bound states like charmonium, only the production of the cc 
pair can be described in pQCD, whereas the formation of the J/i[> bound state which occurs at long distances has to 
be described by phenomenological models. In the framework of non-relativistic QCD (NRQCD) [3j so-called colour 
singlet (CS) and colour octet (CO) states coexist. In the spectroscopy of heavy quark bound states containing only 
one heavy quark, Heavy Quark Effective Theory (HQET) is used to predict the production and decay properties. 
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3. BEAUTY IN PHOTOPRODUCTION 

An important test of pQCD is provided by open heavy flavour pro- 
duction. In Figure [1] the measurements of the beauty cross section 
in photoproduction of both the HI and the ZEUS collaborations are 
shown as a function of the b quark transverse momentum, pr. The 
analyses used different datasets, techniques and decay channels pro- 
viding independent measurements cross-checking each other. The 
different measurements agree well with each other and are in rea- 
sonable agreement with the NLO prediction from FMNPlQ. They 
cover a wide range in p h T giving a consistent picture of b quark pho- 
toproduction. 



4. INELASTIC J/ip PRODUCTION 

At HERA, the charmonium state J/ip is produced predominantly by the boson gluon fusion (BGF) 
process. The theoretical models which are used to compare with the measurements follow either the 
DGLAP evolution or use the kx factorisation approach. After the heavy quark pair (cc) is pro- 
duced at short distances, the formation of the J/tp bound state is described by non-perturbative 
long distance matrix elements (LDME). In the NRQCD models, contributions from both CS and CO 
states are predicted; one aim of the measurements is the extraction of their relative contributions. 
The HI collaboration measured inelastic J /ip production with decays 
to in the photoproduction region using ~ 166pb _1 from the 

2006-2007 data, complementing their previous measurement in the 
DIS region In Figure [2] the differential cross section as a function 
of p T is shown and compared with Monte Carlo predictions and 
the CS calculations at leading order (LO) and at next-to-leading 
order (NLO). The scaled CASCADE Q Monte Carlo, which uses 
the kr factorisation in the colour changing flavour mode, is able to 
reproduce the slope better than the scaled EPJPSI Monte Carlo, 
following the DGLAP evolution. While the CS LO calculation is not 
able to describe the data, there is a good agreement with the NLO 
calculation within the large normalisation uncertainties. As the CS 
provides a generally good description of the data when using the kx 
factorisation or calculations at higher orders, no significant colour octet contribution is required, although there are 
still large normalisation uncertainties on the prediction. 

In order to reduce the effect of the normalisation uncertainty, polarisation measurements are an important testing 
ground for the NRQCD predictions. Normalised quantities can be measured using helicity parameters to characterise 
the decay angular distributions. A measurement of the J/tp helicity [?| has been performed by the ZEUS collaboration 
using the complete HERA statistics (~ 470pb _1 ). The angular distributions of the J/ip — ► decay can be 

parametrised as cx 1 + X(y) cos 2 O + fi(y) sin 20 cos <j> + \v{z) sin 2 6 cos 2 (j) |8|- The following integrated helicity 
formulae are used, depending on the chosen reference frame: 

i- . cx l + A(z)cos 2 ©, 

a a cos Sdz 1 v / ' 

or I + hM. + ±l£i cos 2 tA 

a d<t>dz 3 + 3 Loh 

where 9 is the angle between the \x + vector in the J/ip rest frame and the z axis, and <p is the azimuthal angle in 
the x-y plane of the ji + vector in the J/tp rest frame. 



HERA 





.... i f 


dc/dp?(ep^ebX) 1 




Q 2 <1GeV z , 0.2<y<0.8, |ri b |<2 




EEB NLOQCD(FMNR) 






pg-l/4(lf^ + p3 ] 


7 


H1 99-00 b jet t 


| k T fact. (LZ J2003 set 1 ) '■ 


j T 


H1 99-00 b->ji jet 


I • 


ZEUS 96-00 b-^i jet 




: A 


ZEUS (prel) 05 b-.u jet 




A 


H1 (prel) 06/07 b-»H jet 




o 


ZEUS 96-97 b-»e 




r O 


ZEUS 120 pb 1 b->e 




i □ 


H1 97-00 b->D*(i 




■ 


ZEUS 96-00 b-*D*H 




* 


ZEUS 114 pb" 1 bjb->m» 





5 10 15 20 25 30 



<P?> (GeV) 

Figure 1: Cross sections for beauty production 
as a function of from various decay channels. 




Figure 2: Differential cross section as a function 
of Pt,^ in the ipfp) regime 
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Figure [3] shows the measurement of the helicity parameter v as a 
function of the inelasticity z. In addition to the NRQCD predic- 
tions two different predictions following the k t factorisation approach 
are compared with the measurement. The LO NRQCD predictions 
(BKV) s| do not describe the dependency of this parameter well, 
whereas the prediction including the CO contribution seems to be 
favoured within the large uncertainties. The other two predictions 
(Baranov) [9j], predict a small negative polarisation and are always 
below the measurement. In this scheme where only CS contributions 
have been taken into account two different parametrisations of un- 
integrated gluon distributions have been used. The CS model with 
kx factorisation gives predictions that are more similar to CS+CO 
compared to the CS only prediction. To distinguish between the 
different theoretical contributions the large variation between the 
calculations has to be understood and it would be good if the theo- 
retical predictions could be improved and made available at NLO. 
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Figure 3: Helicity parameter, v, as a function of 
the inelasticity, z 
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5. SPECTROSCOPY 

The large charm production cross section at HERA permits mea- 
surements to be made of excited charm and charm-strange mesons. 
The production of excited charm and charm-strange mesons is ob- 
served by the ZEUS collaboration using the full HERA I dataset 
corresponding to an integrated luminosity of ~ 126pb _1 10]. The 
following decay channels were investigated: 
L>i(2420)° 
L>*(2460)° 

D*+K° s 
-> D*°K+ 

Figure GO shows the M(D* + w a ) and M(D + TT a ) distributions for the 
charm meson candidates reconstructed in the given decay chan- 
nels. The measured masses of the observed mesons have been found 
to be in reasonable agreement with the world average values. In 
addition to the masses also the widths, helicity and the relative 
branching fractions could be extracted. The measured width 
is T(D^) = 53.2±7.2(stat.)±4 g(syst.) MeV which is above the world 
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Figure 4: M(D* + -K a ) and M(D + -K a ) distribu- 
tions for the Di(2420)° and the L>|(2460)° can- 
average value of 20.4 ± 1.7 MeV. A larger S-wave admixture at didates. 
ZEUS with respect to that in measurements with restricted phase space could explain the differences as already 
a small S-wave admixture could have sizeable contributions to the width. The measured value for the he- 
licity parameter of h(D®) = 5.9l^'7(stat.)+i'o(syst.) has to be compared with the prediction of the HQET for a 
pure S-wave (h = 0) and a pure D-wave (h = 3). For the charm-strange meson the measured parameter is 
h(D^i) = — 0.74]to'i7(stat.)]to'o5( s y s t-) which is inconsistent with the prediction for a pure D-wave and more than 
two standard deviations away from the prediction for a pure S-wave. So the measurement suggests a significant con- 
tribution of both D-wave and S-wave amplitudes to the Z? s i(2536) + — > D* + Kg decay. In addition, a search for the 
radially excited charm meson D* (2640) ± has been made, which was reported by DELPHI as a narrow resonance in the 
final state D* ± tt + tt~ at 2637 MeV. In the inspected mass range no signal was observed and according to the expected 



mass and width the following limit at (95% confidence level) was extracted: f(c — > D* + ) • B 



D*' + -*D* + Tt + Tt- 



< 0.4%. 
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6. D*p RESONANCE 

The HI collaboration observed a narrow signal in the D*p chan- 
nel at MJD *p) = 3099 ± 3(stat.) ± 5(syst.) MeV in the HERA I 
dataset ll| . This signal was interpreted as an anti-charm baryon 
with a minimal constituent quark composition of uuddc with a 
relative contribution of ~ 1% in the identified D* me- 

son sample. In contrast to this, no evidence was found by the 
ZEUS pj] collaboration. Using the full HERA dataset HI does 
not confirm the observation. The data of 2004 to 2007 corre- 
spond to an integrated luminosity of 348 pb -1 and thus increase 
the available statistics by a factor of ~ 4. In contrast to the pre- 
vious analysis of the D*p final state no dE/dx requirements were 
used for particle identification. A detailed description of the event 
selection can be found in The mass spectrum was recon- 

structed both for the HERA II dataset as well as for the HERA I 
dataset using the new selection criteria. Figure [5] shows the re- 
constructed invariant mass for the D*p candidates separately for 
the two datasets. In the reanalysed HERA I dataset the previ- 
ously reported signal is reproduced. Under the assumption of the 
width from the HERA I measurement there is no excess in the 
HERA II data visible. The upper limit at a 95 % confidence limit 
was calculated to be ^^.^ ~ 0.1%. 
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Figure 5: Distribution of M(D*p in the HERA I 
(bottom) and HERA II (top) data samples. 
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